Renormalization group approach to anisotropic superconductivity 
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The superconducting instability of the Fermi liquid state is investigated by considering anisotropic 
electron-boson couplings. Both electron-electron interactions and anisotropic electron-boson cou- 
plings are treated with a renormalization-group method that takes into account retardation effects. 
Considering a non-interacting circular Fermi surface, we find analytical solutions for the flow equa- 
tions and derive a set of generalized Eliashberg equations. Electron-boson couplings with differ- 
ent momentum dependences are studied, and we find superconducting instabilities of the metallic 
state with competition between order parameters of different symmetries. Numerical solutions for 
some couplings are given to illustrate the frequency dependence of the vertices at different coupling 
regimes. 

PACS numbers: 74.20-z,74.20.Fg,74.25.Kc 



I. INTRODUCTION 

Many-body interactions have important effects on the 
electronic properties of correlated materials and are sub- 
ject of great attention. The main role played by the 
electron-phonon interaction in the conventional super- 
conductivity has enhanced the interest in electron-boson 
coupling. In the formation of Cooper electron pairs in 
conventional superconductors, the attractive pairing in- 
teraction is mediated by phonons, as established by the 
Bardeen Cooper Schriffer (BCS) theory of superconduc- 
tivity. The Eliashberg theory^'^ provides the appropriate 
equations to obtain the superconducting temperature at 
which pairing occurs, as well as the energy gap created 
in the electronic density of states (DOS). The Eliash- 
berg equations describe an effective electron-electron in- 
teraction due to the exchange of any form of bosons. 
The Eliashberg function a^F(aj), which gives the spec- 
tral electron-phonon density, and the Coulomb pseudo 
potential /i* (Anderson-Morel potential^), which charac- 
terize the pairing interaction, are input parameters of the 
Eliashberg equations. 

Direct evidence that the bosons mediating the attrac- 
tive interaction, in conventional superconductors, are 
phonons was provided by tunneling experiments^'^'^. Fol- 
lowing McMillan and Rowell work^, by inversion of 
tunneling data, a unique value of a^F^uj), as well of 
fj,* , is provided by the structure in the electron tunnel- 
ing current measured as a function of the applied volt- 
age, the dl/dV characteristic. For BCS superconductors, 
the gap created in the DOS when pairing takes place 
is isotropic, of s-wave symmetry. However, anisotropic 
electron-phonon coupling has been reported, later on, in 
some borocarbide material superconductors as YNi2B2C 
and LuNi2C™i^. Furthermore the superconducting en- 
ergy gap of MgB2 was found to vary on the Fermi 
surface (FS) due to the momentum dependence of the 
electron-phonon interaction. The superconducting prop- 



erties of this material, which presents the highest super- 
conducting transition temperature (Tc = 39K) among 
the binary compoundsi^iii, are explained by the fully 
anisotropic Eliashberg theory^il^, by including the mo- 
mentum dependence of the electron-phonon coupling, 
combined with density functional calculations^^. 

The anisotropic Eliashberg theory has also been ap- 
plied to materials with other pairing symmetric a ^"^i^^ such 
as the copper oxides which present a dx^-y^-wave order 
parameter. In this case, besides the phonons, other kind 
of bosons have been considered as, for example, spin 
fluctuations^-. In the copper oxides high- Tc supercon- 
ductors, there is currently not a consensus about the na- 
ture of the pairing interaction. The d-wave symmetry gap 
experimentally measured, seems to favor a purely elec- 
tronic pairing interaction. However a kink, i.e., a change 
of the slope in the quasiparticle energy dispersion, has 
been reported by ARPES experimentsii at an energy 
of about 50meV. This kink could indicate a renormal- 
ization effect similar to that appearing in conventional 
superconductivity. In the cuprates, at this energy scale, 
there are multiple excitations, such as phonons and spin 
fluctuations, which could be responsible for the renor- 
malization. At higher energy, 350 meV, other kink has 
been reported related with bosonic excitations^^. Both, 
phonons-^" and spin fluctuations^^, have been proposed 
but the nature of the bosonic mode remains under de- 
bate. Therefore strong electron-electron correlations as 
well as strong electron-boson coupling have to be consid- 
ered in order to analyze the experimental datap2£. 

Experiments carried out by different techniques on 
hole-doped cuprates, with oxygen isotope substitution 
(i60-i8Q)2j^i22|23^ ji^ye claimed that phonons play an im- 
portant role in cuprates. Recent data have provided ev- 
idence that the electron-phonon interactions are respon- 
sible for the origin of the nodal kink^**, confirming ear- 
lier ARPES results2^ which had indicated an anisotropic 
electron-phonon interaction in cuprates in both the nor- 
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mal and superconducting states. However, the mecha- 
nism of pairing in high temperature superconductors is 
still under discussion. 

We focus here in the study of the superconducting in- 
stability with pairing mediated by anisotropic electron- 
boson coupling. The most general case of anisotropic 
coupling of electrons to bosonic modes is considered, 
with a combination of many symmetry channels. We 
follow the asymptotically exact renormalization group 
(RG) approach developed in refi^. This extension of 
the RG treats electron-electron and electron-phonon in- 
teraction on an equal footing. Analytical and numerical 
solutions of the flow equations for the BCS vertices are 
obtained. Generalized Eliashberg equations and the cor- 
responding MacMillan-Rowell and Allen-Dynes expres- 
sions for Tc are derived from the RG flow equations. The 
Migdal theorerri^l assumed in the Eliashberg theory, is 
understood in terms of a large- expansion in the RG 
approach--. We include both, retardation effects and the 
presence of multiple energy scales in the problem. Re- 
tardation effects in the one-dimensional (ID) Holstein- 
Hubbard model at half-filling have been found to be im- 
portant near the transition between spin density wave 
(SDW) and charge density wave (CDW) states^. The 
classical and quantum aspects of fermion driven lattice 
instabilities have been analyzed using an extension of the 
RG method similar to the one followed in the present 
work^^. The calculation of the gap and the structure of 
the phase diagram of ID tight-binding models, such as 
molecular crystal and Su-Schricffcr-Heeger models, have 
been analyzed at the one-loop levcP-. Furthermore, RG 
methods have been recently used to study the Cooper 
instability in graphene^i^. 

The paper is organized as follows. In section II the 
RG analysis is developed. In Section III the generalized 
Ehashberg equations are discussed . Section IV addresses 
results obtained for electron-boson couplings with differ- 
ent angular dependences. Finally in Section V the work 
is summarized and some conclusions are given. 



a ring of radius kp and width of 2 A around the Fermi en- 
ergy Ep. In the second step, the equations are solved by 
a large- A^ method, with A^ = Ep/ K, N being the number 
of patches at the Fermi surface. Here a 2D square lattice 
is considered, at low fillings, so that the FS has an almost 
circular shape. The interactions are parametrized by the 
on-site Coulomb repulsion uq, the electron-boson cou- 
pling and the Einstein frequency ujp of the bosons. In 
this case the electron-boson constant is A = 2N{0)g^ /lue, 
where A^(0) is the density of states at the Fermi level. 

The RG equations are derived in the path-integral 
representation'^—. In order to obtain a combination of 
many symmetry channels, the most general anisotropic 
electron-boson coupling g{i,j), is considered. The 
electron-boson coupling can be integrated out exactly, 
leading to an electron-electron effective problem with re- 
tarded interactions. The retarded electron-electron in- 
teraction has the form: 

m(4, 3, 2, 1) = u(4, 3, 2, 1) - 2g(l, 3).g(2, 4)D{1 - 3) (1) 



where the phonon propagator is 



(2) 



We adopt the notation q = (i^,q), 1 = (wi,ki), and so 
forth. Spin indices are omitted in this notation. Par- 
ticles 1 and 2 are incoming and scatter into 3 and 4, 
respectively. We will focus on processes involving parti- 
cles with opposite spins. The processes involving parti- 
cles with same spins can be obtained from these due to 
SU(2) symmetry^- . For our almost circular FS, the ef- 
fective retarded electron-electron interaction u gives two 
types of scattering in the RG: the forward scattering, 
with fci = fcs, fc2 = which, as in the case of pure 
electron-electron interaction, does not get renormalized 
but contributes to the electron self-energy, and the scat- 
tering in the Cooper channel. While the forward channel 
does not flow under the RG, the BCS vertex (fci = — fc2, 
/c3 — — /c4) flows at one-loop in the RG equation: 



II. RENORMALIZATION-GROUP ANALYSIS. 

The RG method for interacting electrons^ is adopted 
here to analyze the instabilities of a two-dimensional (2D) 
Landau Fermi liquid. We start from a given geometry of 
the Fermi surface of the non-interacting system and we 
treat the microscopic interactions within one-loop RG, 
which captures the essential physics of the problem. In 
this way, assumption of a predetermined order parameter 
is not needed. Since we are interested in the supercon- 
ducting instability, both electron-electron and electron- 
boson interactions have to be considered. Therefore we 
follow the asymptotically exact RG scheme extended to 
include interacting fermions coupled to bosonic modes^^. 

In the Wilson-like RG theorySS,^ at a flrst step, the 
modes of the momentum space above a cutoff A are inte- 
grated out. The remaining phase space in a 2D system is 



^^(1,3,£): 



duj f'^'^de Aeu{l, uj, 9, l)u{uj, 9, 3, 1) 



27r 



A? 



(3) 



where 1 = {uji,9{) and 3 = (c^s, 9^), and ^ represents the 
RG scale £ = ln(Ao/A), with Aq the initial bandwidth. 
Eq. ([3]) is obtained from the cutoff independent condition 
imposed to the interaction vertex 



m(4,3,2,1) = 0. 



(4) 



The vertices are labeled according to the angle around 
the FS 9i, since the dependence on the radial part of the 
momentum is irrelevant here22. However, the dependence 
on the modulus |q| is important at larger flllings^!^^ 
where, in addition, electron-electron interactions may 
lead to renormalization of the FS shape'^-. In order to 
solve the differential equation ([3]), the pairing potential 
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is decomposed in terms of the irreducible representation 
of the space group of the underlying lattice. For the 
case of a 2D square lattice, the space group is D4. The 
irreducible representations of the D4 space group con- 
tains four one-dimensional (singlets): Ai and A2 for con- 
ventional and unconventional s-wave channels, Bi for 
d2,2_j^2 -channel and B2 for dxy-channel; and one two- 
dimensional (triplet): E, corresponding to p-wave sym- 
metry, with degenerate eigenvalues for the two channels 
Px and py . The corresponding basis functions are chosen 
to be: 





4i{e) - 


Cn 


cos[4(n - 1)61] 


Bi 


€{0) = 


Co 


cos[(4(n - 1) + 2)61] 


B2 


€{o) = 


Co 


sin[(4(n - 1) + 2)61] 


A2 


€{o) = 


Co 


sin[4n0] 


E 


/ €(0) 




Cosin[(2n- 1)6*] 






Co cos[(2n - 1)6*] 



(5) 

with n ranging from 1 to cx), and the normalization fac- 
tors C„ = 1 / \/tt{1 + Sni)- Therefore we can express the 
pairing potential in the more general form 



7 m , Tl 



where 7 labels the representation of the group and 



7 ( 

nin \ 



(6) 



(7) 



where 0^(0) are the corresponding basis functions given 
in Eq. ([5]). The normalized basis functions 0m(0) are 
orthogonal 



27r 



(8) 



and form a complete basis set. By discretizing the fre- 
quency integral into a sum, taking advantage of the or- 
thogonality and completeness of the (/)-basis, and restrict- 
ing ourselves to frequencies below the cutoff, the two- 
dimensional matrix u{i,j,i) expressed in Eq. 16]) can 



be transformed into a four-dimensional tensor 



7^ 



Therefore we can rewrite the flow equation, Eq. 
a tensor equation: 



as 



-[ 



nm\ [3 



M 



N 



where we have imposed an upper limit to the harmonic 
index n < M, and to the frequency sum, —N < 6 < N. 
In Eq. ([9|) we have used 



u 



■7 

nmJ (3 



d9^ d9f3u{a,p,£)^l{9„)(^Z^{00), 



{£) = mLo^j)si 



(10) 



with the kernel 



27r 

d9- 



Therefore Eq. Q can be written in a compact form 
as a {2N + 1)M x (27V + 1)M matrix equation. The flow 
equation have the form 



"df 



-U'^-K^-U'^ 



(12) 



where the matrix indices j are related to the frequency 
a, and harmonic n indices, by a = IP [{j — 1)/Af] + 1 
and n = j — {a — l)M, respectively, where IP denotes 
the integer part. From Eq. (|12p we obtain one vertex 
flow equation for each channel 7. Once we have the ma- 
trix form of the RG flow equations for the vertices, it is 
possible to write the exact solution: 



U^(^) = [1 + U''(0) • Vii)]-^ ■ U'^(O) 



(13) 



where P''(^) = de'K'' {£'). There is an instability of 
the Fermi liquid state when U'^(£c) 00 ioi £ ^ £c- 
This condition is fulfilled when 

det[l + U''(0)-P^(;c)] = 0, (14) 

which is equivalent to solving the eigenvalue equation 

[H-U'^(O) •P''(4)] • =0 . (15) 

The kernel K.'^{£) that appears in the vertex flow equa- 
tions is given by Eq. (|lip and contains self-energy cor- 
rections. The momentum-dependent imaginary part of 
the self-energy E (wjk), with contributions from all the 
components of the electron-boson coupling, is related to 
the quasiparticle-weight by 'Ei{lu,9) = [1 — Ze(Lu,9)]uj. 
The imaginary part of the self-energy is renormalized by 
the bare forward vertices, and the flow equation for the 
quasiparticle-weight Zi{ijj,9) can be integrated to give: 



Zl{uJa,9a) — 1 



dKf, 



d9i3duj(3 



N(Q)uo{a,l3)Z,{u;p,9i3)LO0 

Al + Z^,{LU(3,9f3)Lul ' 



(16) 



where ito(a,/3) is the retarded electron-electron interac- tion 

uo{a,P) ^uo-2g^9^,9p)D{uJa-LO0). (17) 
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The self-energy is angle- dependent and it appears as a 
term in the denominator of K'''(^). Therefore, whereas 
only the 7-component 'K'^{£) of the kernel contribute to 
a given channel 7, all 7-components of the microscopic 
electron-boson coupling contribute to K'' {£) . The results 
of this section will be used in Sec. IV, where a numerical 
analysis of the RG flow equations, as well as the X~Wc 
phase diagram are studied, Wc being the energy scale of 
the instability. 

III. ELIASHBERG EQUATIONS. 

Eliashberg theoryi^ of superconductivity was origi- 
nally formulated to describe phonon-mediated s-wave su- 
perconductors. By including momentum dependence in 
the electron-boson interaction, the anisotropic Eliash- 
berg formalism is obtained. The Eliashberg function 
a^F^Lo), which defines the electron-boson coupling A — 
2 F {Lo)dijj I Lo carries, in the anisotropic formalism, 

momentum dependence. The fully anisotropic Eliashberg 
formalism has been successfully used to study the super- 
conducting properties of MgB2. By solving numerically 
the anisotropic Eliashberg equations, the superconduct- 
ing transition temperature, the momentum-dependent 
superconducting energy gap, and the momentum- 
dependent specific heat have been obtained^'^''^^. Gen- 
eralization of the Eliashberg equations to d-wave super- 



conductivity has as well been done-® and has been used to 
investigate the physics of cuprates^. As it was shown in 
Refi^, the reduced Coulomb repulsion ^* or Anderson- 
Morel potential emerges naturally from the RG equations 
as ^* — uq/(1 -\- £eUo)i where £e = ln(Ao/ci;£;). The 
value of Tc, and the analytical forms of the McMillan^ 
and AUen-Dynes^i^ expressions are obtained from the 
RG flow equations. 

Here the instability conditions obtained in the pre- 
ceding section for the BCS vertex are mapped into the 
Eliashberg equations. The breakdown of the Fermi liq- 
uid state by the SC instability, due to the retarded 
electron-electron interactions, occurs at a finite tempera- 
ture which turns out to be the same critical temperature 
Tc obtained from the Eliashberg theory. To investigate 
this instability, we solve the RG flow equations at finite 
temperature. Since the SC instability is approached by 
decreasing the temperature, in the equations such as (jlSp . 
where integration over A is involved, this can be extended 
to zero. In this formalism, considering a finite value of 
the temperature T, the quasiparticle weight Zi{io, 6) also 
presents a T dependence. Taken the appropriate integra- 
tion limits in Eq. (|16p we obtain a similar expression to 
Eq. (jlSp . which gives the condition for the appearance 
of the instability and allows us to obtain the expression 
for Tc- The integrals over frequencies can be written 
as Matsubara sums (/ duj0/{2Tr) TcJ^p) i'^ equations 
(fT5)) and leading to: 



Z{ujn,e,Tc) = 1 + 7rrc^fe(w„ - w,„). (18) 



These two expressions are the set of generalized Eliash- 
berg equations at Tc, where 

fe^ (iu^-ujp) = -( ^7io(a, /?) (19) 

with Ti^ {9(3,9) = E^i'/'?(^/3)'^?(^)- The uo term that 
appears in the definition of Tg^ via ug, does not con- 
tribute to Z because it does not have any dependence on 
frequency, and therefore the Matsubara sum in Eq. (fTS]) 
vanishes for this term. The mq term in Eg g, does however 
contribute to the Eq. involving v'^{LLimO)- There is 
one such generalized Eliashberg equation for each chan- 
nel 7, which depend on the 7-component of the initial 
coupling through T^g g, as defined in Eq. pop . However, 
each of these equations, corresponding to a particular 
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channel 7, also depends on the quasiparticle weight Z 
which is renormalized by all the components of the initial 
coupling through Tg as defined in Eq. (fT9)) . 

From the generalized Eliashberg equations, one can ob- 
tain the value of Tc in terms of the microscopic parame- 
ters. For the most general case, analytical expressions for 
Tc are quite involved. We consider first the special case 
in which g{9a,9i3) is separable, so that it has the form 
g{()a,9p) — f{9a)f{()i3), and set ito = for now. Then 
Tg g, will also be separable: 

Tl g,{uJa - t^/3) = 'i^gi'^a " U!p)T2,{uJa - t^/s), (20) 

where 

T2iLo^-L,p)= f f\9^)^2NlS))D{Lo^ - ujp)7f {9^,0). 



5 



On the other hand, it is always possible to write: 

fg{uJa = ^ J " ^^/3)f - UJp), (21) 

where Tg{uja — w^) = J^j^ei^a ~ ^l^)- general 
case, since we know that g{9a,dp) — g{9j3,9a), we can 
use the decomposition g{9a,9f5) = J2i fii(^a)fi{d(3), and 
therefore: 

i 

T0{uJa-COp) = 77- / T ie {UJ a - OJ f3)r le '{UJ a -Lup) (22) 

The set of equations derived here reduce to those given 
by earlier treatment of Eliashberg theory for anisotropic 
electron-boson interaction by Daams and Carbotte^^, if 
we consider a single separable channel, 7 = 1 and n = 1. 

If for a given channel 7, the gap is dominated by 
one component s such that the SC gap can be written 
as A^(0) = Ao(j)']{9), then the corresponding McMillan 
expression^ for Tc becomes: 



T|? w uje exp 



fi*6. 



7I '^7 



(23) 



where d-yi restricts the contribution of /x* to the channel 
7=1 {uq only has first harmonic of s-wave component), 
Z-i=/,0^(0)/[l + A(0)] and 



\{9^,9p)c^l{9^)7f(9p,9)c^1{9) 
l + \{9) 



,(24) 



with \{9a,9f3) = 2N{0)g^{9a,9p)/ujE, and X{9) = 
{2tt)^^ Jg, X{9,9'). The corresponding Allen-Dynes ex- 
pression is 



(25) 



where 

p 



A(0„, 0^)02 (0„)77^ (0^5, 0)0^(0). (26) 



The effective electron-phonon parameter [Eq. ([M)) ] 
that appears in the McMillan and Allen-Dynes equations 
for a particular channel 7 does include contributions from 
all components of the bare coupling A(0q, 9^) through 
and X{9). 



IV. NUMERICAL RESULTS 

The method used allows for an easy numerical calcu- 
lation of the critical temperature (as well as the zero 
temperature gap Aq) from the microscopic parameters 
Mo, X{9a,9f:j) and uje- Eq. gives the RG evolution 
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FIG. 1: Plots of the {2N + 1)M x (2iV + 1)M matrix Vie) 
at different RG scales £ for the s-channel (7 = 1), for the 
electron-boson coupling g{0,6') — gocos{9 — d'). Here, for 
the first harmonic in the expansion (M = 1), the number of 
frequency divisions is 41 (N — 20) and Aq — 100, ue ~ 10, 
uo — 0.1. The three panels on the top correspond to A = 0.4 
(weak coupling) and the three on the bottom correspond to 
A = 4 (strong coupling). 



of all the vertices and can be evaluated at different RG 
steps £. To find the instability, one can simply calculate 
the quantity det[l -|- U'*'(0) ■ P^(£)] and see when it ap- 
proaches zero. Here we study electron-boson couplings 
with different angular dependences and strengths, going 
from the weak to the strong coupling regimes. We la- 
bel 9 the angle associated to k = ki = — k2 in polar 
coordinates, and 9' the angle of k' = ka = — k4. In gen- 
eral we define g{9, 9') = gaf{9, 9'), where 50 is a position 
independent constant and the function f{9,9') contains 
the full angular dependence of the electron-boson ma- 
trix elements. We consider two different cases: a) The 
electron-boson coupling only depends on the difference 
between the angles of the electrons, b) The coupling de- 
pends on the exact position of each of the the electrons 
on the Fermi surface. 



A. Solution of the RG flow equations for 

g{0,9') = 5ocos(6i - e') 

For this kind of angular dependence, the strength of 
the electron-boson coupling does not depend on the spe- 
cific position of the particles on the FS, but only on 
their relative orientation. By solving Eq. (fT4)) we ob- 
tain the evolution of 'U"'{£) with £. For the coupling 
g{9,9') — gocos{9 — 9'), only three channels contribute, 
corresponding to s-, dx^-y^- and c?a;y-symmetry, respec- 
tively. Furthermore, the contributions of the two chan- 
nels with d-symmetry have the same magnitudes, the 
channels are degenerate, and we will generically refer 
to them as d-wave channel. In Fig. [T] the evolution of 



6 



■ ■ ■ ■ -20 W ■ ■ ■ ■ -20 ' 



-20 -10 10 20 -20 -10 10 20 -20 -10 10 20 




-20 -10 10 20 - 20 -10 10 20 - 20 -10 10 20 



FIG. 2: Same as Fig. [T]but for the d-channel (7 = 2, 3). The 
three panels on the top are for weak coupling (A = 0.4) and 
the ones on the bottom are for the strong coupling regime 



(A = 4). 



Xi'^ii') with ^ for the s-channcl is depicted. Each panel 
represents the (2iV + 1)M x (27V + 1)M metric U^(.^) 
at a given RG step t. The matrix elements correspond- 
ing to small frequencies, around = 0, are at the cen- 
ter of the panels. The three top panels of Fig. [1] show 
the RG flow evolution for a case in the weak coupling 
regime (A = 0.4). The first panel on the left shows the 
initial condition, for which £ = 0. The third one, on 
the right hand side, shows the vertices at £ = £c, right 
before the instability. We see that the most divergent 
couplings are those with frequencies below the Einstein 
frequency, Iwaj, jw/jj < ue- In this case, simple two-step 
RG can be applied. The three panels on the bottom of 
Fig. [1] represent the evolution of the couplings for the 
strong electron-boson interaction regime (A = 4). In this 
case, the most divergent couplings (at I = ^^ are those 
with frequencies below a given energy scale Wc > i-^e- 
As expected, these results are similar to that obtained 
in Ref."*^ where an isotropic electron-boson coupling was 
considered. Furthermore, this can be seen as a proof of 
the validity of our method, which matches the results of 
the isotropic limit when only the s-channel is considered. 

The situation is quite different for the d-channel, as it 
can be seen in Fig. [51 Again, the panels on the top corre- 
sponds to the weak coupling regime. The initial condition 
is the same as for the s-channel, with the u'^{u!, ~lo, £ = 0) 
being the largest couplings. However, close to the insta- 
bility, the couplings u''(wq., oj/j, ^c) that first diverge cor- 
respond to a region of frequencies centered around lo ~ Q 
and with a star-like shape. This frequency dependence 
of m'^ at the instability point illustrates the importance 
of retardation effects in the RG flow. To have a physical 
understanding of this speciflc structure is, however, not 
trivial, and we do not speculate here about the possible 
physical consequences of such dependence. On the other 
hand, the bottom panels of Fig. [5] represents the flow. 



in the strong coupling regime, from the initial condition 
^ = 0, to the critical RG step Ic- In this case, the most 
divergent couplings at the instability point are, as in the 
s-wave case, those with frequencies below Wc- 

In Fig. [3] we show the density plots of the matrix 
U*(£c) at the RG scale ic, where the instability ap- 
pears, for different values of the electron-boson coupling 
strength. At this point we can see more clearly how 
the energy scale separating high and low energy physics 
moves from the Einstein frequency w^;, in the weak cou- 
pling regime, to the critical cutoff Wc, in the strong cou- 
pling regime. The scale Wc can be associated with the 
T = superconducting gap Aq, or with the critical tem- 
perature Tc of the SC phase in the finite temperature 
formalism2S'43, pjg [3] also illustrates the breakdown of 
the two-step RG. In this approximation, the vertex is cho- 
sen to be just the electron-electron part for frequencies 
above the Einstein frequency, and to have some constant 
contribution from the boson modes for frequencies below 
loet^- This approximation works well in the weak cou- 
pling limit, where the most divergent couplings are those 
with frequencies below lue- But for large A (A ^ 1) this 
behavior breaks down and the scale for the divergent cen- 
tral region is of order Wc > uje- 

Similarly, the couplings U'^(£c) in the d-channel are 
shown in Fig. [51 from the weak to the strong electron- 
boson coupling regimes. We obtain the star-like struc- 
ture at weak coupling (A = 1), evolving towards a more 
circular form in the strong coupling (A = 6). More ana- 
lytical study is required for a full understanding of these 
patterns. 

Fig. [51 shows the energy scale Wc, where the SC insta- 
bility occurs, as a function of the strength of electron- 
boson coupling A, at fixed loe- The two channels (s 
and d) that contribute to the electron-boson coupling 
g{9,9') = gocos{9 — 0'), are represented. The behav- 
ior of Wc in the weak coupling limit follows the McMil- 
lan exponential behavior of Eq. (|23p . For large A, strong 
coupling regime, the Wc follows the Allen-Dynes law, Eq. 
(|55)) . In the inset of Fig. \5\we show the phase diagram 
of the system. At finite temperature, the Fermi liquid 
phase breaks down towards d-wave SC in the weak to in- 
termediate coupling range, and towards s-wave SC in the 
intermediate to strong coupling regime. The larger the 
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FIG. 3: Couphng matrix U''(£c) at the instability point for the 
s-channel, from weak to strong coupling, for uq — 0.1. The 
corresponding electron-boson coupling is g{9, 6') — go cos{6 — 
9'). Panels correspond, from left to right, to A = 1,2.5,6. 
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on-site electronic repulsion wq, the larger the electron- 
boson coupling crossover between s- and d-wave super- 
conductivity. This is expected since an isotropic uq only 
suppresses the contribution to the s-wave channel, with- 
out affecting the other channels. Therefore for uo = 0, 
superconductivity only occurs in the s-channel. Notice 
that, because the Kohn-Luttinger effect is not included in 
this one-loop RG calculation, there is no superconduct- 
ing instability when only electron-electron repulsion is 
considered. In addition, this effect should shift the tran- 
sition between d-wave and s-wave superconductivity to- 
ward some different value of the electron-boson coupling, 
A. The qualitative behavior is, however, well captured 
by our approximation, as shown in the phase diagram of 

Fig.m 



B. Solution of the RG flow equations for 

5(6', 61') = 3ocos(6i/2) cos(6l72) 




In the following we consider an electron-boson 
coupling with the angular dependence g{9,9') = 
gocos{9/2)cos{9' /2). Unhke the previous case, this is a 
simple example for which the boson mode couples differ- 
ently to electrons in different parts of the Fermi surface. 
That is, the coupling depends on the specific position of 
each of the electrons in the Fermi surface. As before, we 
obtain the -^-evolution of XJ'^{i) by solving Eq. (fTl)) for 
this coupling. The two contributing channels in this case 
are of s- and p-symmetry. The evolution of \Jp{£) for 
the p-channel, at weak electron-boson coupling, shows a 
star-shaped structure, as it was found in the preceding 
section for the d-channel. The frequency dependence of 
the coupling matrix close to the instability is very similar 
to that represented in the top right panel of Fig. [21 

In Fig. [HI the density plots of ^{Ic) are shown for 
Mo = 0.1 and different values of A, at the critical RG step 
ic- They present the same qualitative behavior as the 
corresponding s-channel results of the previous subsec- 
tion (see Fig. 2]) . For the p-channel, the density plots of 
the coupling matrix U^(^c) at the instability point £c, are 
depicted in Fig. [T] The plot for A = 1.0 is very similar to 
that obtained for the d-channcl, in the case of a coupling 
of the form g oc cos{9 — 9'). This can be observed by com- 
paring the left hand side graphs of Fig. |3]and[71 However, 
the situation is different in the strong coupling limit. The 



FIG. 5: (Color online) Critical cutoff Wc versus A for the s- 
channel (black squares) and the d-channel (red circles). The 
value of the parameters used are Ao = 100, loe = 10 and 
Uo = 2.0. The inset represents the phase diagram. 

'jl HI iHII 'jBUl 

-10 ^^^^1 ^^^^1 -10 ^^^H J^^^H -10 ^^^m 

-20 -10 10 20 -20 -10 10 20 -20 -10 10 2( 

FIG. 6: Coupling matrix U^(^c) at the instability point 
for the s-channel, from weak to strong coupling, for uo — 
0.1. The corresponding electron-boson coupling is g{0,0') = 
go cos(S/2) cos(&'/2). Panels correspond, from left to right, to 
A = 1,2.5,6. 



corresponding plot for the d-channel shows that the most 
diverging matrix elements lie in a circular shaped region 
around a; = (see the right panel of Fig. 0]). However 
in the present case, where g oc cos(0/2) cos(0'/2), the di- 
verging matrix elements for the p-channel are confined in 
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FIG. 4: Same as Fig. [3] and for the same parameter values, 
but for the d-channel, V^iic). 



FIG. 7: Same as Fig. [6] and for the same parameter values, 
but for the p-channel, U''(^c). 
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a square region around w = 0, as shown in Fig. [T] This 
different behavior in the strong couphng region needs a 
better understanding. 

Finally, we mention that it is possible to build a phase 
diagram similar to that of Fig. [5l but for a coupling of 
the form g cx cos(0/2) cos(0'/2). In this case, instead 
of having a low-coupling d-wave SC region in the phase 
diagram, we obtain a zone with p-wave pairing SC. For 
a value of mq = 0.5, the crossover between p- and s-wave 
superconductivity occurs at A ~ 2.5. For a small value 
of the on-site Coulomb interaction, mq = 0.1, it is found 
that the Fermi liquid state is unstable towards p-wave 
SC in the range < A < 0.7, and towards s-wave SC for 
larger values of the couphng, A > 0.7. 

V. CONCLUSIONS 

In this work, we have investigated the pairing insta- 
bilities of the FL state which appear when the electron- 
boson coupling overcomes the effective repulsive electron- 
electron interaction. The most general case of anisotropic 
coupling of electrons to bosons has been considered. The 
anisotropic boson exchange couplings are treated by an 
RG approach. Calculations are made for a 2D square 
lattice at low fillings. Under this assumption, the system 
presents an almost circular FS, which allows for an ana- 
lytical solution of the RG flow equations. Furthermore, 
the numerical evaluation of the flows provides informa- 
tion of the change of the frequency dependence of the ver- 
tices. The flow equations for the BCS vertex are decom- 
posed into contributions coming from different symmetry 
channels, with different angular momentum dependence. 
By varying the strength of the interaction, from the weak 
to the strong coupling regime, the evolution of the cou- 
plings for the different symmetry channels is obtained. 
Channels of s- , p-, and d-symmetry have been investi- 
gated. 

We have considered here simple functional forms of the 
electron-boson coupling, such as g{9, 9') — go coa{9 — 9'), 
where the momentum of the exchanged boson only de- 
pends on the angle difference between the two electrons 
involved, and g{9,9') — go cos{9/2) cos{9' /2), where the 
coupling depends explicitly on the position of each of the 
electrons on the Fermi surface. It is found that, even 
with these simple angular dependences, the anisotropic 
electron-boson couplings induce new non-trivial physics. 
As far as we know, the inclusion of angular depen- 
dence of the electron-phonon coupling has not been ad- 
dressed before in previous RG calculations. Although 
academic at first view, this problem makes contact with 
the physics which appears in some new anisotropic mate- 



rials as those described in the Introduction. In these ma- 
terials, metal-transition borocarbides, boronitrides, mag- 
nesium diboride, or cuprates, standard BCS theory can- 
not explain many of their properties. A consistent de- 
scription of some physical behaviors is however achieved 
if anisotropy of the electron-boson coupling is included 
in the Eliashberg theory^i^i^^. In high-temperature su- 
perconductors, there is a considerable amount of data 
that point to interplay between electronic and atomic 
degrees of freedom^. The different behavior shown by 
the quasiparticles in both nodal and antinodal regions of 
the Brillouin zone, has suggested, among other explana- 
tions, an anisotropic electron-phonon coupling, along side 
other many-body effects, in order to understand the pair- 
ing mechanism^'*®. Coupling to the half-breathing mode 
in the nodal region and to the buckling mode in the 
antinodal direction have been propose as an interpreta- 
tion to the renormalization effects seen in ARPES results 
in both, the normal and superconducting phases^S.. How- 
ever, more work is needed to get insight in the complex 
interplay between electron-boson interactions and elec- 
tronic correlation in unconventional superconductivity. 

In summary, we have studied the superconducting in- 
stability of the Fermi liquid phase, considering electron- 
boson couplings with different angular dependences and 
strengths. SC order parameters of s, p and d symme- 
try have been obtained, depending on the anisotropy of 
the coupling and the strength of the interaction. The in- 
vestigation of the frequency dependence of the couplings 
of different symmetries could be interesting in order to 
analyze and understand the complex behavior revealed 
by the experiments. Therefore, numerical evaluation of 
the RG flow equations are instructive in determining the 
range of important frequencies in different regimes and 
for different symmetry channels. Furthermore, at the in- 
stability point and for finite temperatures, the RG equa- 
tions give the solution of the generalized Eliashberg equa- 
tions at Tc and consequently, McMillan and Allen-Dynes 
expressions have been obtained. 
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